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ABSTRACT 


The three-dimensional flow in the region where a turbulent boundary layer encounters an 
appendage is characterized by the formation of a junction root vortex. The formation of this 
vortex can be attributed to the adverse pressure gradient upstream of the appendage. One 
major area of concern is the temporal behavior of the vortex. Fluctuations in the size and/or 
location of the vortex can cause degradation of downstream components or excitation of 
local structures. Recent studies in wind tunnels have identified a bi-modal characteristic in 
the structure of the vortex. 


An experimental investigation of the time-dependent behavior of the junction root vortex was 
conducted using MIT’s variable pressure water tunnel and laser doppler anemometry system. 
Data was collected at various points forward of a modified NACA 0020 section appendage 
mounted on a splitter plate and analyzed using off-line DFT analysis. The data was gathered 
at tunnel speeds of 10 and 20 ft/sec using an analog-to-digital acquisition system. Velocity 
probability distributions using a fully digital acquisition system were also plotted to identify 
the bi-modal structure. 


The spectra generated from the DFT analysis were inconclusive. No distinct frequencies or 
pattern of frequencies attributable to the flow could be identified. Fluctuations within the 
spectra made their accuracy questionable. The existence of the bi-modal structure was con- 
firmed by the velocity probability distributions. 
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1 INTRODUCTION 


1.1 BACKGROUND 


The two-dimensional turbulent boundary layer along a body has been studied exten- 
sively and is well understood. When this boundary layer encounters an appendage, how- 
ever, the flow in the region of the junction becomes extremely complex. This 
three-dimensional flow is characterized primarily by the formation of a "horseshoe" or 
junction root vortex at the base of the appendage. There are many practical examples of 
this flow including the keel of a sailboat, a submarine’s sail and control surfaces, an air- 


plane’s wings and stabilizers, and even a telephone pole in a windstorm. 


The formation of this vortex secondary flow can be attributed to the adverse pres- 
sure gradient upstream of the appendage. The vorticity in the incoming turbulent bound- 
ary layer’ is concentrated or coalesced by the gradient forming the vortex. The flow 
causes the vortex to wrap around the appendage turning the cross-flow vorticity of the 
incoming boundary layer into streamwise vorticity. The vortex produces a three- 


dimensional separation region in the vicinity of the appendage. 


Dickinson [5] divided the flow around an appendage into three conceptual regions. 
The inflow region is the region upstream of the appendage where the presence of the 
appendage is not felt. The second region is the appendage region. This is the region 
which includes the adverse pressure gradients and the generation of the secondary flows 
associated with the appendage. This region is further divided into the leading edge subre- 
gion, the alongside subregion, and the trailing edge subregion. The final region is the 


wake. 


1 A junction vortex can also form when the incoming boundary is laminar. 
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A number of experimental studies have been conducted or are currently in progress 
to characterize the junction vortex. Most of the early work concentrated on investigating 
the flow for various appendage geometries. Love [9] conducted some preliminary 
research which showed that the curvature of the leading edge affected the strength of the 
generated vortex. Since then, Belik [2] and Baker [1] have investigated the flow using 
cylindrical appendages; while Shabaka and Bradshaw [21], Oguz [13], Kubendran [8], 
and McMahon et al. [10],[11] conducted studies using infinite-chord appendages. To 
study the flow around more realistic hydrodynamic shapes, Dickinson [4],[5] of the 
David Taylor Research Center (DTRC) and Pierce et al. [14],[15],[16] of Virginia Poly- 
technic Institute (VPI) have used streamlined foils in their experimental work. A more 
complete literature review on the topic of the junction vortex can be found in Sung [22]. 
This document not only references experimental studies, but also publications on analytic 


and numerical work. 


The majority of experimental work on the junction vortex flow has been for the pur- 
pose of obtaining data to support the development of analytic or numerical prediction 
schemes. The data has been collected using various experimental techniques including 
flow visualization, hot-wire/film anemometry, laser doppler anemometry, and pressure 
taps with most of the investigations having been conducted in low-turbulence wind tun- 
nels. Recently, a great deal of attention has been focused on reducing the strength of the 
generated root vortex. In marine applications, the presence of this vortex can degrade the 
performance of downstream control surfaces and propellers. The use of fillets and other 
leading edge modifications is currently under study by various institutions. Pierce, Fran- 
gistas, and Nelson [17] have published some preliminary results from the studies being 


conducted at VPI. 
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Another area of concern is the temporal behavior of the junction vortex. Fluctua- 
tions in the size and/or location of the vortex can cause additional degradation of down- 
stream components, or excitation of local structures resulting in vibrational and acoustical 
problems. Rood et al. [18],[19] conducted some of the first studies on the tume-dependent 
nature of the junction vortex flow. He reported relatively low frequency and narrow 
bandwidth turbulent fluctuations in the three-dimensional boundary layer downstream of 
a streamlined protuberance with a 1.5:1 elliptical nose, similar to the one used in this 
study. In a meeting on appendage design at DIRC in 1988, Dr. R.L. Simpson [3] of VPI 
reported detecting a bi-modal characteristic in the structure of the vortex. Histograms of 
streamwise velocity fluctuations measured close to the line of shear* exhibit two distinct 
peaks. Pierce et al. [17] have shown that the spectra within the vortex flow is similar to 
that in the upstream boundary layer, suggesting that the fluctuations in the large scale vor- 


tex structure originate in the boundary layer. 


1.2 RESEARCH GOALS 


This experimental investigation had two purposes. The primary goal was to gather 
raw data on the junction vortex flow and process it to obtain information on the time 
dependent nature of the flow through spectral analysis. Unlike previous studies which 
were conducted using wind tunnels, this investigation was conducted in a circulating 
water tunnel. An attempt to verify the bi-modal structure identified by Simpson and 


quantify its frequency of fluctuation was a principal element of this objective. 


2 The line of shear delineates a region of high shear near the foil from an area of lower shear outside. The line 
of shear falls inside the separation line forward of the appendage. 
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The other purpose of the investigation was to determine the feasibility of using 
velocity data from the laser doppler system in MIT’s variable pressure water tunnel facil- 
ity to conduct spectral analysis. A key element of this goal was to incorporate a new 


all-digital acquisition and flow analysis system into the laboratory facility. 
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2 EXPERIMENTAL APPARATUS 


2.1 VARIABLE PRESSURE WATER TUNNEL 


All testing and data acquisition was conducted using MIT’s variable pressure water 
tunnel [7]. This facility is designed to test and evaluate the dynamic performance of pro- 
pellers and two and three dimensional hydrofoils. The test section is a 54 inch long, 20 


inch square closed-jet section. The maximum test section speed is 33 ft/sec. 


The area directly upstream of the test section contains the upstream turning vanes, a 
circular to square transition section, a settling section, and a contraction. The transition 
section has a contraction ratio of 1.46:1, while the ratio for the contraction section is 
4.92:1. The overall contraction ratio for the facility is 7.17:1. The settling section con- 
tains a honeycomb consisting of one inch diameter acrylic tubes, each one being one foot 
long, and three wire screens. The settling section is designed to reduce the turbulence in 
the test section. Measurements of the free stream turbulence using the laser doppler velo- 
cimeter show a turbulence intensity of less than 0.5 percent, accounting for phase noise’ 


in the velocity measurements. 


The test section is immediately followed by a diffuser, downstream turning vanes, 
and a square to circular transition section. The closed circuit flow in the water tunnel is 


driven by a 30.25 inch diameter four-bladed impeller. 


2.2 BODY/APPENDAGE JUNCTION MODEL 


The model used in the study of the junction vortex consisted of a modified NACA 


0020 foil section mounted on a flat splitter plate. The foil and plate were both 


3 Phase noise is related to the bandwidth of the doppler signal spectrum and can be related to an equivalent 
turbulence intensity. 
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constructed using anodized aluminum. The fou section was produced on a numerically 


controlled mulling machine. 


Figures 1 and 2° show the model layout, test section mounting arrangement, and 
gu y £ g 


experimental coordinate system. 
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Figure 1. 
Test Section Layout 
(Front View) 


4 All dimensions in Figures 1 and 2 are in inches. 
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TEST SECTION WINDOW 
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Figure 2. 
Test Section Layout 
(Side View) 


The shape of the foil section was chosen to conform to several previous investiga- 
tions (see Section 1.1). The standard NACA 0020 section was modified by replacing the 
nose with a 1.5:1 elliptical arc. This bluff leading edge produces a stronger adverse 
pressure gradient which promotes the formation of the junction vortex. Equations 2.1 and 
2.2 were used to generate the foul section offsets. 


Nose Section: 


; = 18) (2.1) 
DEG Cen oe 
12° ° ) kes 
Tau Section: 
yie = 0.2969 - Vx/c -0.126 - x/c — 0.3516 - (x/cy (a) 


+0.2843 - (x/c)Y —0.1015 - (x/cy" 
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The two shapes were joined at their respective maximum widths (2.4 inches) 
resulting in a chord length of 10.3 inches. The span of the foil was set at 15 inches to 


allowing a 2.4 inch gap between the foil and top of the test section. 


HALF-THICKNESS (IN ) 





CHORD LENGTH (IN.) 


Figure 3. 
Modified NACA 0020 Foil Section 
Tnp wires of circular cross-section were positioned just aft of the leading edge of 

the splitter plate and 0.4 inches upstream of the maximum thickness of the foul. The trip 
wire on the splitter plate invoked the development of a turbulent boundary layer along the 
plate forward of the foul. As recommended by Devenport and Simpson [3], the boundary 
layer on the foul was tripped to avoid any unsteadiness of the flow which might have 
resulted from natural transition. A trip wire diameter of 0.047 (3/64) inches was used in 
both cases making them fully effective down to a tunnel speed of 2.26 ft/sec as stated by 


Schlichting [20]. 


erie) Ne 





2.3 LASER DOPPLER VELOCIMETER SYSTEM 


The velocity data used for this study was gathered using the horizontal component 
of a three-component laser doppler velocimeter. A three watt Argon-Ion laser serves as 
the coherent light source for the measurement system. The beam of light from the laser is 
split into three components, green (A=5.14 nm), blue (A=488.0 nm), and violet (A=476.5 
nm), by an optical color separator. The blue component is used for horizontal measure- 


ments. 


Each component of the laser is further split into two beams in a dual beam system. 
By use of optical lenses and murors, the two beams intersect at a fixed angle inside the 
tunnel test section forming an ellipsoidal measurement volume. The two intersecting 
beams of light produce a fringe pattern within this measurement volume perpendicular to 
the plane formed by the two beams. A particle in the fluid moving through this fringe 
pattern will scatter light. The intensity of the scattered light will be modulated at a fre- 
quency determined by the fringe spacing and velocity of the particle perpendicular to the 
fringe pattern. This frequency is called the doppler frequency and is given by Equation 


BS. 


2-u-sink (2.3) 
pe <= 


Jp = doppler frequency 
u = component of velocity perpendicular to the fringe pattern 
k = one-half of the intersection angle of the two beams 
X = wavelength of light 


The scattered light from particles crossing the measurement volume 1s collected by 
a set of receiving optics and detected with a photomultiplier tube which converts the 


modulated light signal to a voltage signal with the same doppler frequency. The signal is 
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then processed by a data acquisition system to obtain the measured component of the par- 
ticle velocity. The receiving optics can be placed at any location where the measuring 
volume can be effectively viewed. In this study, the same optics were used to focus the 


outgoing light and collect the scattered light; this is called an on-axis backscatter system. 


A Stationary fringe pattern cannot differentiate between equal positive and negative 
flow velocities, nor can it convey a zero flow condition. To overcome this shortcoming, a 
frequency shift is applied to one of the two crossing beams. This causes the fringe pattern 
to slowly move upstream or downstream depending on the magnitude and direction of the 
shift. The moving fnnge pattern will produce a doppler frequency even in a zero flow 
condition and will provide different frequencies for equal positive and negative flow 


velocities. 


The laser doppler velocimeter used in this study was mounted on a three-axes tra- 
verse table. The motion of the table was controlled remotely from a personal computer 
using feedback control. The positional accuracy of the traverse table, as stated in the 


operations manual, 1s 2.5 um. 


Additional information on laser doppler anemometry theory and techniques can be 


found in Durst, Melling, and Whitelaw [6]. 


Table I. 
Laser System Characteristics 


188.0 am 
Diameter of 63 um 
Measurement Volume 

Length of 736 um 
Measurement Volume 


Number of Sa 







Page 17 





3 DATA ACQUISITION 


3.1 ANALOG-TO-DIGITAL ACQUISITION 


The data for the spectra in Appendix B was collected using an analog-to-digital 


acquisition system. Figure 4 depicts this acquisition system. 
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Figure 4. 
Analog Data Acquisition System 





As shown, the modulated voltage signal from the photomultiplier tube was fed to an 
analog tracker. The tracker is a signal processor specifically designed to convert the 
doppler frequency to an analog voltage, and extract the signal from the background noise. 


Since the signal from the photomultiplier tube is discontinuous and varies in amplitude, 
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the tracker serves as a combination of tracking filter, frequency to voltage converter, and 
sample and hold circuit. The validity of the signal is confirmed before making it avail- 


able for output as a voltage directly proportional to flow velocity. 


The output voltage from the tracker was sent through a signal conditioner. The only 
purpose of this device was to apply a 1kHz low pass filter to the voltage signal, which 
was necessary to prevent aliasing during the off-line spectral analysis (see Section 4.1). 
The filtered voltage was then converted to digital information by using one channel of a 
16-channel analog-to-digital converter board installed in the laboratory’s personal com- 


puter. 


The operation of the converter board was controlled via software directly from the 
computer keyboard, permitting selection of sampling frequency and number of data 
points. The digitized voltage versus time data was converted to an output file and placed 
on a floppy disk for later analysis. For each collection of samples, a graphic display of 
voltage versus time was available on the computer screen. Figure 5 is an example of a 


typical display. Note the large amount of scatter in the tracker output. 
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Voltage 





Time (sec) 


Figure 5. 
Typical Digitized Data Display 


The digital-to-analog converter and data analysis software used to generate the 
spectra in Appendix B was designed for operation on continuous signals. The tracker’s 
output signal, however, appears as "points" or "lines" due to the data validation circuity. 
This simply means that the last reading is held until a new reading is validated. The out- 
put then is basically a series of point measurements rather than a continuous output. To 
meet the continuity requirements for the analysis, it was necessary, therefore, to insure 
that the data rate from the tracker sufficiently exceeded the data sampling frequency. In 
all cases, the tracker data rate exceeded the 2000 Hz sampling frequency by a factor of 


five. 
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3.2 FULL DIGITAL ACQUISITION 


The data for the velocity probability distributions in Appendix C was collected 


using a fully digital acquisition system. Figure 6 depicts this acquisition system. 
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Figure 6. 
Digital Data Acquisition System 





The primary difference between the data acquisition system in Section 3.1 and that 
in Figure 6 1s incorporation of a digital signal processor or Intelligent Flow Analyzer 
(IFA). One exceptional feature of the IFA is that it makes only "good" measurements or 
it makes no measurements at all. The IFA determines the doppler frequency by measur- 


ing the duration of eight contiguous cycles of the doppler signal. This elapsed time is 


5 The Intelligent Flow Analyzer (IFA) is a product of Thermal Systems, Inc. (TSI). 
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then processed to obtain the frequency. To extract the true doppler signal from 
background noise, the IFA uses an autocorrelation of zero-crossings technique. The 


result is highly accurate data. 


The IFA was not used to gather the data for the spectral analysis for two reasons: 
(1) It was not available during the initial measurements. (2) The processor data rate was 
insufficient to meet the 2000 Hz sampling requirement. In many cases, the maximum 
obtainable data rate was lower by a factor of ten than that required by the spectral analy- 
sis (see Section 4.1). Non-optimal tunnel seeding® and the processors stringent validation 


algorithm were the probable causes of this unfortunate limitation. 


The IFA was directly controlled from the laboratory personal computer. The real- 
time display software installed on the computer not only controlled the IFA functions, but 


also produced the velocity probability distributions found in Appendix C. 


6 Seeding is the process of adding particles to the flow to improve the fluids ability to scatter light. 
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4 DATA ANALYSIS 


4.1 SPECTRAL ANALYSIS FUNDAMENTALS 


The purpose of spectral analysis is to determine the frequency composition of a 
continuous random function of time x(t). Since most experimental measurements of ran- 
dom functions are now conducted digitally, the spectral analysis is performed on a dis- 
crete series of numbers, x,(t), obtained by sampling the continuous function x(t). The 
introduction of the fast Fourier transform (FFT) has made the calculation of spectral 


estimates directly from the discrete time series extremely quick and very accurate. 


Any non-periodic function of time can be represented by a Fourier integral of the 


form 
°° °° 4.1 
x(t) = { A(w) cos(wt )dw +f B(w)sin(wt)dw Sey 
w = angular frequency (rad/sec) 
where 
l f) 
A(w)=— | x(t)cos(wt)dt 
ae J = (4.2) 


oo 


B(w)= | x(t) sin(wt )dt 


—oo 


For engineering purposes, Equation 4.1 is true in classical Fourier analysis if 


. 4.3 
{ Ix(t) | dt <0 ne 


—o0 


Defining X(u') as 


X(w) =A(w) —iB(w) (4.4) 
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Equations 4.2 can be combined to give 


~ 4.5 
X(w) =< | x(nje"at ae 


Equation 4.5 is the formal definition of X(w) and 1s called the Fourier transform of x(t). 


The condition given by Equation 4.3 cannot be satisfied by an unending stationary’ 
process. Therefore the classical Founer theory cannot be applied to most experimental 
data. This difficulty is overcome by analyzing the autocorrelation function R,(t). The 
autocorrelation of a function x(t) is defined as the average of the product x(t)x(t+t), which 
for a Stationary process 1s only a function of the time difference t. The autocorrelation 
function can be adjusted to meet the condition of Equation 4.3 allowing it to be Fourier 


transformed. The transform of R,(T) 1s called the spectral density or the spectrum of x(t). 


Most experimental data, as stated above, 1s a discrete series of samples x,(t), r= 
0,1,2,...,.(N-1), where t=rA and A=T/N. T is defined as the total sampling time and N is 
the total number of samples. In this situation, the integral 1n Equation 4.5 can be replaced 


by the summation 


1 A 4.6 
X.=— econ Wia=(0.102...,(N = 1) oe 


Equation 4.6 is the formal definition of the discrete Founer transform (DFT), where & 1s 
the harmonic number that is related to a discrete angular frequency w = 27%k/T. The fast 


Founer transform is a computer algonthm for calculating DFT’s. 


7 A process is said to be stationary if the probability distribution of a collection of samples is independent of 
absolute time. 
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The advantage of the DFT is that it can be used to generate the spectrum of x(t) 
directly from the discrete time senes x,(t) without first generating the autocorrelation 
function. The theory supporting this can be found in Newland [12]. The DFT produces 
spectral components at the frequencies defined by the harmonic numbers k = 
0,1,2,...,(V-1). Values of the spectrum for k > (N-1) can be shown to simply repeat them- 
selves. Additionally, it can be shown that the spectrum is symmetrical about the zero 
frequency position. The unique and valid portion of the spectrum exists only between 
w = 0 and w =7/A. The upper limit on the frequency is called the Nyquist frequency’ w, 
and is the maximum frequency that can be detected from data sampled at a time spacing 
A. If frequencies greater than w, exist in the original signal, distortion of the spectrum 
called aliasing can occur. The original signal must be filtered prior to processing to 


remove these frequencies to avoid this distortion. 


The accuracy of a spectral measurement depends on the effective bandwidth of the 
measurement and the total sampling time T. The ratio of the standard deviation 6 to the 


mean value m of the spectral density 1s given by Equation 4.7. 


1 (4.7) 


B.T 


Cael 
—_ 
— 


x |aQ 


For most practical methods of calculating the spectrum using DFT’s, the effective 
bandwidth B, is approximately 1/f Hz. This means that 6/m = 1, or the magnitude of the 
standard deviation of the spectrum is equal to the mean. To improve this accuracy, 
several adjacent spectral estimates can be averaged increasing the effective bandwidth. 


Equation 4.7 becomes 


8 The Nyquist frequency, when converted to Hz, is f, = 1/2A. 
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m_ N 2n +1 


(2n+1) = number of adjacent estimates 
The larger the number of adjacent estimates that are averaged, the greater the accuracy, 
but not without a price. The increase in effective bandwidth also lowers the frequency 


resolution. 


Another method of improving the accuracy of a spectral estimate is the application 
of a tapered data window. This function modifies the input data record by weighting the 
middle of the digital information more heavily than the two ends. The overall effect is a 
smoothing of the input data with a corresponding increase in the sharpness of the spectral 
window. The increase in sharpness comes at the expense of a slightly larger bandwidth 
and lower frequency resolution. Figure 7 shows one type of a tapered data window called 


a Hamming window. 
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Figure 7. 
Hamming Window 
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The following steps, as given by Newland [12], summarize the procedure for devel- 


oping spectra from digital input data: 


1. Estimate the frequency range of interest. 

2. Choose an appropriate sampling interval A to insure that the Nyquist fre- 
quency exceeds the frequency of interest. Filter the input data to remove fre- 
quencies higher than the Nyquist frequency. 

3. Decide on the required accuracy 6/m. 

4. Estimate the required frequency resolution and specify the maximum effec- 
tive bandwidth B.,. 

5. Calculate the total sampling time T from Equation 4.7. 

6. Determine the number of data points needed from N = T/A. The total num- 
ber of data points for each record must be equal to a power of two to apply 
FFT analysis. 

7. Determine the number of adjacent spectral estimates that must be averaged 
from Equation 4.8. 

8. Conduct the Fourier analysis using data taper windows if desired. 

9. Average the required number of adjacent spectral estimates to obtain the 


smooth spectrum. 


4.2 SPECTRUM DEVELOPMENT 


The spectra in Appendix B were generated using the theory and procedure outlined 


in Section 4.1. Table II lists the particular parameters used in the Fourier analysis. 
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Hable i: 
Fourier Analysis Parameters 


Frequency Range 0-500 Hz 
Nyquist Frequency f, 1000 Hz 
Sampling Frequency 2000 Hz 
Sampling Interval A 0.0005 sec 


Number of Sample Points NV 4096 





The data for the spectra was gathered by surveying the area just forward of the foul, 
the incoming boundary layer, and the free stream. Measurements were made at tunnel 
speeds of 10 and 20 ft/sec. This equates to Reynolds numbers of 8.7x10° and 1.7x10° 
based on the chord length of the foul. If Reynolds numbers are calculated using the 
momentum thickness of the theoretical boundary layer, as in previous studies, the results 
are 5330 and 9280 respectively. Attempts were also made to collect data at a tunnel 


speed of 25 ft/sec, but tracker lock-on could not be maintained. 


The coordinate system shown in Figures 1 and 2 was used as a reference for the 
data collection. The origin of the coordinate system was at the junction of the foul and the 
splitter plate as indicated in Figure 2. The location of each measurement point was given 
in millimeters to allow use of the traverse table software. The area forward of the sail 
was surveyed on the plane of symmetry using a two dimensional grid of six horizontal 
points and three vertical points. The size of the grid was set to incorporate the significant 


portion of the junction vortex. The boundary layer was sampled at a position 255 mm 
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upstream of the foil, the same location used to measure the boundary layer profiles given 
in Appendix A. The free stream data was collected in the test section with only the split- 


ter plate installed. 


The data collected by the analog-to-digital acquisition system was placed on floppy 
disks for off-line analysis. The first step of the analysis was to put each data file through 
a digital 4th-order Butterworth bandpass filter. For most of the spectra in Appendix B, 
the bandpass filter was set at 10-1000 Hz. The upper limit was set to match the setting of 
the analog low-pass filter of the signal conditioner in the acquisition system. Once again, 
the value of this cutoff was determined from the Nyquist frequency. The lower limit was 
chosen to remove low frequency fluctuations in the voltage signal caused by fluctuations 


in the tracker data rate. Figures 8 and 9 show the effect of the bandpass filtering. 
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Figure 8. 
Unfiltered Data 
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Figure 9. 
Data after 10-1000 Hz Bandpass Filter 


Figure 8 was produced using the same data as Figure 5, except the data has been time 
averaged to reduce the clutter of the display. Note the fluctuation in the voltage starting 
at approximately 0.5 seconds. This was caused by a drop in the tracker data rate. As can 
be seen in Figure 9, the bandpass filter removed the fluctuation providing smooth data for 
the spectral analysis. Also note that the mean value or 0 Hz component of the voltage 
was removed. This was of no consequence since the spectral analysis software would 


have removed this component prior to generating the spectrum. 


The filtered data was then used to generate the power spectral density (PSD) for 
each measurement point. In all cases, the spectra were normalized to a maximum value 


of 1.0. Because of this normalization, some of the data files were filtered using a 50-250 
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Hz bandpass filter and then analyzed to enhance the spectral peaks in the 100 to 200 Hz 
range. These spectra can also be seen in Appendix B. A Hamming tapered data window 
was applied to all data files during the analysis to enhance and sharpen the spectral peaks. 


A typical output of the spectral analysis is shown in Figure 10. 





0.5 
+ 
> 
oe | 
i | 
2s i 
| 
o4 
+ 
+ 
: | 
0.3 - 


PSD 





Frequency (KHz) 


Figure 10. 
Power Spectrum using Hamming Window 


The final step necessary to produce the smooth spectra in Appendix B was to per- 
form the averages of the adjacent spectral estimates. Five adjacent values were averaged 
producing the accuracy given in Table II. Figure 11 depicts a final smooth output 


spectrum. 
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Figure 11. 
Smooth Power Spectrum 
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5 RESULTS 


Appendix A contains the measured profiles of the approach boundary layers for tunnel 
speeds of 10, 20, and 25 ft/sec. As stated before, these profiles were measured at a position 
255 mm upstream of the leading edge of the foil on the plane of symmetry. The purpose of 
these measurements was to insure that a turbulent boundary layer was properly developing 
along the splitter plate. The flow was tripped just downstream of the leading edge of the 


plate to induce this development (see Section 2.2). 


The curves shown in Figures 12, 13, and 14 are the profiles predicted theoretically 
assuming a 1/7th power law for the velocity distribution and the measured thickness of the 
boundary layer. Each experimental data point was determined by averaging the values of 
1024 samples taken over a 200 sec tume period. The data shows a great deal of scatter, espe- 
cially at 10 ft/sec, much more than would be expected. The higher the flow rate, the closer 
the data approaches the predicted profile. In all cases, the measured data generally follows 
the predicted velocity profile; therefore, the boundary layers were accepted as satisfactory for 


the remainder of the study. 


The first group of spectra in Appendix B are those for the free stream at 10, 20, and 25 
ft/sec. At least two measurements were made and analyzed at each speed to determine the 
reproducibility of the spectra from the digitized analog data. These spectra also serve as a 


baseline to which the spectra in the boundary layer and junction region can be compared. 


The consistency of the free stream spectra at 10 and 20 ft/sec is limited. In both cases, 
the approximate size and number of frequency peaks are similar between measurements, but 
the location or frequencies of the peaks appear shifted. The magnitude and direction of the 
shift varies depending on the frequency, possibly indicating that the source of these various 


frequencies is not common to all. The majority of the peaks do not shift by more than 50 Hz, 
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but there are some exceptions. If the first spectrum from 10 ft/sec is neglected, the remaining 
two are almost exactly alike with the exception of the frequency shifts, especially in the 500 


Hz to 1000 Hz region. 


The collection of spectra generated using the data from the boundary layer and junction 
region are extremely disappointing. They appear to contain very little, if any, quantitative 
information. In view of the inconsistency of the free stream spectra, no particular frequen- 
cies can be isolated with any degree of accuracy. All attempts to correlate a spectrum within 
the junction region with other spectra in the region, the boundary layer, or the free stream 
prove unsuccessful. No distinct frequencies or frequency pattern can be identified and 


directly attributed to the flow within the region. 


A small amount of qualitative information may be extracted from the spectra despite 
their limitations. Near the plate (Z = 1mm), at a tunnel speed of 10 ft/sec, the spectrum for 
X = 5mm indicates a very insignificant distribution of energy in the frequency range of inter- 
est (0-500 Hz). This spectral energy distribution changes significantly, however, as the posi- 
tion of measurement point moves away from the leading edge of the foil along the plane of 
symmetry. The spectrum from data taken at X = 10mm shows a large increase in all spectral 
components producing numerous frequency spikes. Based on the predictions of previous stu- 
dies and knowing the location of the flow’s stagnation point from mean velocity measure- 
ments, this position is close to the core of the junction vortex. As the measurement point 
continues to move away from the leading edge, the trend appears to reverse resulting in lower 


average spectral levels as the spectra approach the levels found in the boundary layer. 


It is difficult to confirm the above trend in spectral energy levels from the data gathered 
at 20 ft/sec. In general, the levels measured in the junction region at this tunnel speed are 


significantly lower than those recorded at 10 ft/sec. 
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The velocity probability distributions in Appendix C are much more conclusive in pro- 
viding useful information on the nature of the junction flow. These distributions were 
created using the new digital flow analyzers in conjunction with the laser doppler system. 
The low data rate, which made these analyzers impractical for the spectra development, was 
not a factor in this portion of the study. The IFA’s inherent accuracy, however, was crucial 
for this part of the investigation. The distributions in Appendix C clearly show the bi-modal 


characteristic 1n the structure of the junction root vortex. 


The velocity probability distribution of the free stream at 10 ft/sec is exactly what 
would be expected. The shape is Gaussian with a very small standard deviation. The slight 
unsymmetrical nature of the curve about the mean is caused by the input filter settings of the 
IFA. The turbulence intensity calculated from this distribution is less than 0.5 percent when 
adjusted for the phase noise of the laser signal. The distribution in the boundary layer is also 
a Gaussian curve with a larger standard deviation caused by the increased turbulence in this 


region. 


The velocity probability distributions within the junction region are unlike those of the 
free stream and boundary layer. In all cases, the distributions appear to be the sum of two 
Gaussian distributions displaying two distinct peaks. The relative strength of the two peaks, 
as expected, is dependent on the distance away from the leading edge of the foil. By follow- 
ing the progression of figures in Appendix C, it is obvious that the larger peak shifts from a 
negative velocity behind the stagnation point to a positive velocity as the stagnation point is 
crossed. The presence of two peaks is indicative of two distinct flow velocities existing over 


tume at the measurement point confirming a bi-modal structure. 
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The velocity probability distribution results obtained at a tunnel speed of 20 ft/sec are 
the same as those acquired at 10 ft/sec. The bi-modal characteristic of the flow structure is 
even more pronounced at this tunnel speed. The two distinct peaks are clearly visible at X = 


5mm and X = 10mm. 
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6 CONCLUSIONS AND RECOMMENDATIONS 


The conclusions that can be drawn from this experimental study of the junction root 
vortex are not as extensive as originally expected. The results of the velocity probabulity dis- 
tributions clearly confirm the presence of the bi-modal characteristic within the structure of 
the junction vortex as stated by Simpson. The spectral studies, however, provide no clear 
information about the flow except to indicate an increase in the spectral energy distribution in 
the vicinity of the vortex. The frequency of fluctuation of the bi-modal structure or any other 
distinctive frequency pattern caused by the flow cannot be identified from the spectra 


obtained in this investigation. 


The ability to use data from MIT’s laser doppler anemometry system to obtain spectral 
information is still questionable and requires further study. The accuracy of the data gath- 
ered using the IFA’s could significantly improve the quality of the spectral output if suffi- 


cient data rates could be obtained. 


The following is a list of recommended procedures for further investigation. Some of 
these procedures relate directly to the investigation of the junction vortex, while others 


broadly apply to the use of the laser doppler system for tume-frequency studies: 


1. Conduct on-line frequency analysis of the laser doppler signal using a spec- 
trum analyzer and the analog output of a signal processor. The sample-and-hold 
nature of both the tracker and IFA must be taken into account. A sufficient data 
rate, depending on frequencies of interest, must be achievable to assume a contin- 
uous input signal to the analyzer. 


2. Improve the data rates obtainable with the IFA to allow the collection of 


9 Although the primary output of an IFA is digital information, an analog output is available from a connection 
on the back panel of the processor. 
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accurate digital information for off-line spectral analysis within a useful fre- 
quency range. 

3. If sufficient data rates can be obtained from the IFA’s on both the horizontal 
channel and fiber optic (third component) channel, attempt to cross-correlate 
fluctuations in the junction region with fluctuations in the free stream or bound- 
ary layer. 

4. Improve the accuracy of the spectra by averaging a larger number of adjacent 
spectral estimates. A 6/m of ten percent can be obtained by averaging 100 adja- 
cent values. 

5. Average a large number of spectra generated from sets of data taken at the 
Same point over time. This may help to reduce the fluctuations in the spectra, as 
indicated in the free stream cases, allowing only significant frequency compo- 


nents to emerge. 
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APPENDIX A. 
BOUNDARY LAYER PROFILES 
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Figure 12. 
Boundary Layer Profile (10 ft/sec) 
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Figure 13. 
Boundary Layer Profile (20 ft/sec) 
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Figure 14. 
Boundary Layer Profile (25 ft/sec) 
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APPENDIX B. 
POWER SPECTRAL DENSITIES 
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Figure 15. 
Free Stream, 10 ft/sec 
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Figure 16. 
Free Stream, 10 ft/sec 
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Figure 17. 


Free Stream, 10 ft/sec 
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Figure 18. 
Free Stream, 10 ft/sec 
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Figure 19, 
Free Stream, 10 ft/sec 
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Figure 20. 
Free Stream, 10 ft/sec 
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Figure 21. 
Free Stream, 20 ft/sec 
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Figure 22. 
Free Stream, 20 ft/sec 
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igure 23: 
Free Stream, 20 ft/sec 
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Figure 24. 
Free Stream, 20 ft/sec 
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Figure 25. 


X=5mm, Z=1mm, 10 ft/sec 
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Figure 26. 
X=5mm, Z=1mm, 10 ft/sec 
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Figure 27. 
X=5mm, Z=I1mn, 10 ft/sec 
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Figure 28. 
X=10mm, Z=1mn, 10 ft/sec 
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Figure 29. 
X=10mm, Z=1mn, 10 ft/sec 
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Figure 30. 
X=10mm, Z=1mn, 10 ft/sec 
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Figure 31. 
X=15mm, Z=1mn, 10 ft/sec 
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Figure 32. 
X=15mm, Z=1mm, 10 ft/sec 
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Figure 33. 
X=15mm, Z=1mm, 10 ft/sec 
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Figure 34. 
X=25mm, Z=1mn, 10 ft/sec 
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Figure 35. 
X=25mm, Z=1mm, 10 ft/sec 
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Figure 36. 
X=25mm, Z=1mm, 10 ft/sec 
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Figure 37. 
X=5mm, Z=8mm, 10 ft/sec 
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Figure 38. 
X=5mm, Z=8mm, 10 ft/sec 
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Figure 39, 
X=5mm, Z=8mn, 10 ft/sec 
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Figure 40. 
X=15mm, Z=8mn, 10 ft/sec 
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Figure 41. 
X=15mm, Z=8mn, 10 ft/sec 
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Figure 42. 
X=15mm, Z=8mm, 10 ft/sec 
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Figure 44, 
=25mm, Z=8 


Z=8mm, 10 ft/sec 
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Figure 45. 
X=25mm, Z=8mm, 10 ft/sec 
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Figure 46. 
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Figure 47. 
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X=5mm, Z=2.8mm, 10 ft/sec 
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Figure 48. 
X=5mm, Z=2.8mm, 10 ft/sec 
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Figure 49. 
X=15mm, Z=2.8mm, 10 ft/sec 
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Figure 50. 
X=15mm, Z=2.8mm, 10 ft/sec 
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Figure 51. 
X=15mm, Z=2.8mm, 10 ft/sec 
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Figure 52. 
X=25mm, Z=2.8mm, 10 ft/sec 
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Figure 53. 
X=25mm, Z=2.8mm, 10 ft/sec 
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Figure 54. 
X=25mm, Z=2.8mm, 10 ft/sec 
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Figure 55. 
Boundary Layer, 10 ft/sec 
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Figure 56. 
Boundary Layer, 10 ft/sec 
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Figure 57. 
Boundary Layer, 10 ft/sec 
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Figure 58. 
X=5mm, Z=1mm, 20 ft/sec 
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Figure 59. 
X=5mm, Z=l1mm, 20 ft/sec 
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Figure 60. 
X=10mm, Z=1mm, 20 ft/sec 
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Figure 61. 
X=10mm, Z=1mm, 20 ft/sec 
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Figure 62. 
X=15mm, Z=1mn, 20 ft/sec 
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Figure 63. 
X=15mm, Z=1mn, 20 ft/sec 
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Figure 64. 


X=25mm, Z=Imn, 20 ft/sec 


Figure 65. 


X=25mm, Z=1mmn, 20 ft/sec 
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Figure 66. 
X=5mm, Z=7mn, 20 ft/sec 
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Figure 67. 
X=5mm, Z=7mn, 20 ft/sec 
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Figure 68. 
X=15mm, Z=7mm, 20 ft/sec 
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Figure 69. 
X=15mm, Z=7mm, 20 ft/sec 
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Figure 70. 
X=25mm, Z=7mm, 20 ft/sec 
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Figure 71. 
X=25mm, Z=7mm, 20 ft/sec 
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Figure 72. 


X=5mm, Z=2.1mm, 20 ft/sec 
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Figure 73. 
X=5mm, Z=2.1mm, 20 ft/sec 
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Figure 74. 
X=15mm, Z=2.1mm, 20 ft/sec 
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Figure 75. 
X=15mm, Z=2.1mm, 20 ft/sec 
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Figure 76. 
X=25mm, Z=2.1mm, 20 ft/sec 
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Figure 77. 
X=25mm, Z=2.1mm, 20 ft/sec 
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Figure 78. 
Boundary Layer, 20 ft/sec 
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Figure 79. 
Boundary Layer, 20 ft/sec 
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APPENDIX C. 
VELOCITY PROBABILITY DISTRIBUTIONS 
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Figure 80. 
A beduahra 10 ft/sec 


The PSD is Gau with a small standard deviation 
n-symmetry 1s due Sarees settings of the signal processor. 
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Figure 81. 
Boundary Layer, 10 ft/sec 
The larger standard deviation is due to turbulence in the boundary layer. 
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Figure 8&2. 
X=5mm, Z=l1mm, 10 ft/sec 


Negative velocity is dominant near the foil as expected. 
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Figure 83. 
=10mm, Z=1mn, 10 ft/sec 
The double wee or bimodal structure Starts to become evident. 
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Figure 84. 
X=12.5mm, Z=I1mn, 10 ft/sec 
uble peak is evident. 
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Figure 85. 
X=15mm, Z=1mn, 10 ft/sec 
The dominant velocity shifts to positive flow. 
Bimodal structure is still evident. 
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Figure 86. 
X=20mm, Z=Ilmn, 10 ft/sec 
The dominant velocity is now strongly positive. 
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Figure 87. 
A ete 20 ft/sec 
Gaussian, similar to 10 ft/sec. 
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Figure 88. 
Boundary Layer, 20 ft/sec 


Turbulence in the boundary layer once again broadens the distnbution. 
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Figure 89. 
X=1.5mm, Z=1mn, 20 ft/sec 
Double peak is already starting to become evident very close to foul. 
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Figure 90. 
X=5mm, Z=l1mm, 20 ft/sec 
The bimodal structure is clearly evident with a dominant negative velocity. 
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Figure 91. 
X=10mm, Z=1mm, 20 ft/sec 
The bimodal structure ts strongly inherent at 20 ft/sec. 


Page 95 





4, O8E 





—J 2 Preic soo rio 


6.29 


0.8 


4, 000 





4, Be 


WE 
VELOCITY (M/$) 


Figure 92. 
X=15mm, Z=1mm, 20 ft/sec | 
The dominant velocity shifts to positive. 
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